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ABSTRACT 

By adopting the differential age method, we utilize selected 17832 luminous red galaxies (LRGs) from 
Sloan Digital Sky Survey Data Release Seven (SDSS DR7) covering redshift 0-0.4 based on 'Carson 
& Nichol sample' to measure Hubble parameters. With single stellar population (SSP) models, we 
derive optimal age information of our selected sample. From the decreasing age-redshift relation, we 
obtain the new observational H(z) data (OHD) points which are H (z) — 69. 0± 19.6 km s" 1 Mpc -1 at 
z = 0.07, H{z) = 68.6 ± 26.2 km s" 1 Mpc" 1 at z = 0.12, H(z)=72.9 ± 29.6 km s^ 1 Mpc" 1 at z = 0.2 
and H(z)— 88.8 ± 36.6 km s _1 Mpc" 1 at z = 0.28, respectively. Combined with other 21 available 
OHD data points, a good performance of constraint on ACDM model is presented. 
Subject headings: Cosmology:cosmological parameters -Cosmology: observations -Galaxies: evolution 



1. INTRODUCTION 

The expansion of the Universe can be quantitatively 
studied using the results from a variety of cosmological 
observations, for example the mapping of t he cosmic mi- 
crowave background (CMB ) anisotropies (|Spergel et al.l 
l2007t iKomatsu et al.ll2010l) . the measurement of baryon 
acoustic oscillation ( BAO) peaks (|Eisenstein et alJfeOOSl ; 
IPercival et al.l I2010[). the linear growth of large-scale 
structures (LSS) l|Wang fc Tegmarkl 120041), strong gray - 
itational lensing ([Yang fc Chen! 120091: ICao et al.ll2012h . 
and measurements of 'standard candles' such as the 
redshift-distance re lation s hip of type la supernovae 



(SNIa) (IRiess et al.l U99l ILi et al 



2011 ) and gamma - 
ray bursts (GRBs) (jGhirlanda et al.H2004tiLi et al.ll2008l ). 
Among the various observations is the determination of 
the Hubble parameter H(z) which is directly related to 
the expansion history of the Universe by its definition: 
H(z) = a/a, where a denotes the cosmic scale factor 
and d is its rate of change with respect to the cosmic 
time. The method based on the observational H(z) data 
(OHD) has been used to test cosmological models (e.g. , 
Yi fc Zhandl2007t ISamushia fc Ratrall2006t IZhang et al 



2007T 
2010; 



Zhan g fc Zhuli2 008: Carval ho et aUl2008t IFeng fc li 
Chen fc Ratrall201 1: Xu 2012). Besides parameters 



constraints, OHD c an also be use d as an auxiliary model 
selection criterion (|Li e t al. 2009|). 

In practice, the Hubble parameter H(z) is usually mea- 
sured as a function of the redshift z, with a(t)/a(to) = 
1/(1 + z): 

f > 1 dz 

H ( z ) = —^— — - (!) 



1 



dt 



According to equation ((T|), H (z) has been measured 
through the differential m e thod, which was first put for- 
ward by Uimenez fc Loe b (2002). Derivative of redshift 
with respect to cosmic time, dz/dt has a direct deter- 
m ination on H(z ) . Thi s me thod has been demo nstrated 
in Uimenez et al.1 (j2003T ) and iSimon et al.l (|2005[ ) who ex- 
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tended the determination of H (z) to 8 more redshift bins 
up to z 1.8. Th is dataset was brought up-to-date by 
IStern et al.1 (f2010h . 

However, it may be difficult to determine the accu- 
rate age of a galaxy considering stars in a galaxy born 
continuously, and a young st ellar population m ay dom- 
inate their emission spectra ([Zhang et al.l 120101 ) . Fortu- 
nately, it has been found in great number that luminous 
red galaxies (LRGs) which have photometric properties 
consi stent with an old, pass ively evolving stellar popula- 
tion ([Roseboom et al. 2006) co uld be selected as 'cosmic 
chronometers' ([Crawford et al.l 120101 ) . 

Different from other methods of obtaining the 
age history of galax ies by co-adding spectra (e.g., 
ICarson fc~ Nichol 2010), we employ ULySS, an available 
package on-line to fit a full-length spectrum. Through 
single stellar population (SSP) fitting, age information 
could be obtained which will be discussed in detail in 
Section [31 From the age-redshift relation, OHD points 
could be deduced finally. 

There are already 21 OHD points obtained from both 
differential age method and baryon acoustic oscillation 
(BAO) method (pvforesco et al|[20ll IStern et al.| [2010l: 
ISimon et~a l. 2005; Gaztana ga et al.ll2009D . Currently the 
number of OHD points are still scarce compared with 
SNIa luminosity distance data. Therefore the potential 
power of O HD in con s traini ng cosmological parameter is 
explored in lMa et al.l ([201 ID in detail. It has been found 
that as many as 64 further independent determination 
of H (z) are needed to match the parameter constraining 
power of SN la. In Section [51 we combine our 4 new data 
points with those 21 OHD points to constrain cosmolog- 
ical parameter. 

This paper is organized as follows. We briefly give 
our LRGs selection sample in Section [5] In Section [31 
we describe explicitly how we gain the age information 
from the LRG spectra using ULySS and in Section 2] we 
present our method to get the OHD using the galaxy 
ages. A cosmology constraint using all available OHD 
including our 4 new ones is given in Section [5j Finally, 
in last section, we discuss the limitation and implication 
of our results. 
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Figure 1. Redshift distribution of 17832 LRGs. 

2. SOURCE SELECTION OF LRG SAMPLE 

SDSS DR7 (lAbazaiian et al.l 12009ft includes 11,663 
deg 2 of imaging data, with most of the 2000 deg 2 in- 
crement over the previous data release lying in regions 
of low Galactic latitude. The catalog contains five-band 
photometry for 357 million distinct objects. The sur- 
vey also includes repeat photometry on a 120° long, 2.5° 
wide stripe along the celestial equator in the Southern 
Galactic Cap, with some regions covered by as many as 
90 individual imaging runs. 

The LRG selection criterions of SDSS lEiscnstci n et al.l 
(2001) are on the basis of color and magnitude to yield a 
sample of luminous, intrinsically red galaxies. However, 
the sample selected from SDSS according to Eisenstein's 
selection criterion are not very homogeneous. Our sam- 
ple has been extracted fr om Carson fc Nichol's sample 
(Carson & Nichol sample; iCarson fe~N ichol 2010) with a 
further restraint on signal to noise ratio (SNR). As SNR 
has an impact on our next fitting step, it is necessary 
to demand SNR of R-band to be greater than 10. The 
17832 selected quiescent, luminous red galaxies basing on 
'Carson & Nichol sample' method from SDSS DR7 cover 
the redshift range from 0-0.4, and the distribution of all 
LRGs is shown in Figure [TJ 

Compared with other spectra selection work, 'Carson 
& Nichol sample' calibrates, for the first time, the SDSS 
spectra on Lick/Image Dissector Sca nner (IDS) system 
(|Burtein et al.1 119841 lWorthev et al.l I1994D . Their gen- 
eral selected criterion steps are as followings. First step 
is to get L RGs from Catalo g Arch ive Server (CAS) as 
outlined in lEisenstein et al.l ([2001D . Further, restrain- 
ing the spectrum using standard emission lines such as 
Ha, H/3 and OIII 5007 to increase the number of truly 
quiescent galaxies in the previous sample. To fully pic- 
ture physical properties such as velocity dispersion and 
absolute luminosity of these LRGs, through correcting 
velocity dispersions for aperture effects and performing 
K + e corrections to the magnitude, four subsamples are 



produced with different absolute magnitude and veloc- 
ity dispersion. For a explicit descriptio n of this method, 
please refer to ICarson fe Nichol (|2010D . 

3. AGE- REDSHIFT RELATION 

We proceed to describe our way of obtaining age in- 
formation of LRGs. In this paper we employ ULySS 
(jKoleva et al.l l2009al ) to analyz e the selected spectra. 
This package is available at |http://ulyss. univ-ryonl.fr/} 
In ULySS, an observed spectrum is fitted witfi a model 
expres sed as a linear com bination of non-linear compo- 
nents (jKoleva et al.ll2009al ) . The full spectrum is utilized 
in fittings for determination of the s tar formation an d 
metal enrichment history of galaxies (jWu et al.l l2011af) . 
This package minimizes x 2 between an observed spec- 
trum and a model spectru m. This method has already 
been successfully tested in IWu et all (|2011bf) . The reli- 
ability and robustness of using ULySS to the study of 
the history of stellar population has been verified (e.g., 
i Koleva et aHl2009bt iKoleva et alJl2009d : iBouchard et alJ 
I201CA iMakarova et alJl2010f >. 

3.1. Matching resolution 

The first step in fitting is to match resolutions be- 
tween the observed spectra and the model we use. To 
match the resolution of our LRG spectra from SDSS 
DR7 to the models, we transform the model provided by 
ULySS by injecting relative line spread function (LSF) 
between our spectrum and the model. When deter- 
mining the LSF, we make use of the available veloc- 
ity dispersion templat e stars as a standard star from 
http://www.sdss.org/dr7/algorithms/veldisp.html. We 
make a five times linear interpolation in wavelength and 
convolve with the model. In this way, a new matching 
model is generated, and all this p rocess could be acco m- 
plished through ULySS function (jKoleva et al.ll2009al ). 

3.2. Single stellar population result 

It is essential to study the stellar population of galax- 
ies if we want to reconstruct the star formation history 
(SFH) of galaxies. We fit the spectrum with a single 
SSP for the following reasons. Firstly, LRG is believed 
to be drawn from a same parent population with the 
most percent of their stars are formed from a single- 
burst. Secondly, SSP-equivalent properties correspond 
to the 'luminosity-weighted' average over the distribu- 
tions. Thus by a single SSP fitting, a general view of a 
galaxy could be gotten (jDu et alll2010l) . We carry out an 
initial stu dy of age-dating us i ng th e stellar population 
models of iBruzual fc Chariot! ((2001 (BC03) to synthe- 
size spectra, and this model is also provided by ULySS. 

BC03 is a model for computing the spectral evolution 
of stellar populations at ages between 1 x 10 5 yr and 
2x 10 10 yr at a resolution of 3 A across the whole wave- 
length range from 3200 A to 9500 A f or a wide range of 
metallicities (jBruzual fc Charlotll2003|) . The accuracy of 
this model has been proved and it has been widely used 
for getting information of spectra. But when interpreting 
the galaxy spectra, we find both of age and metallicity 
guess value tend to affect the fitting result. It is diffi- 
cult to remo ve such influen ce especially at the low resolv- 
ing spectra pu et aW20lH . The fitting of ULySS starts 
from a point (initial value) in the parameter space(age 
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Figure 2. Redshift-metallicity relationship from \ 2 map (40x40). 
Open circles correspond to a subsamplc of 164 LRGs. It shows that 
all galaxies metallicity are consistent and fluctuate between 0.1 to 
0.2 dex approximately. 

and metallicity), therefore we stress the importance of 
that point as the returning fitting value will be initial 
value dependent and the minimization will be local. In 
order to identify and understand the age-metallicity de- 
generacy and local minimizations, we analyze that age- 
metallicity degeneracy and construct \ 2 map which is a 
function provided by ULySS. The \ 2 map is a visualiza- 
tion of parameters space, and the principle of this map 
is to give a grid of nodes in a 2-D projection of the pa- 
ramet er space. \ 2 map ret urns the minimization of each 
node ([Koleva et al.ll2009al ). 

A subsample of 164 LRG spectra from redshift 0.03 
to 0.179 is dealt with by building their x 2 maps. The 
metallicity results of this sample are shown in Figure O 
From Figure [5J it has been discovered that the metal- 
licity of LRGs are consistent with the assumption that 
selected galaxies are though t to be have similar metal- 
licities ([Jimenez et alJl2~003[ ) and it floats between about 
0.1 dex and 0.2 dex (see Figure [2]). 

ULySS could provide initial value vectors as initial 
value, that is to say, if giving points intensive enough, it 
is possible to find the global minimization and thus break 
the degeneracy. We set our metallicity range guess points 
from 0.1-0.2 dex to every galaxy in our sample, while age 
guess point span from 5000 Myr-11200 Myr. 

For each galaxy, groups of age and metallicity values 
can be gained as results. We ignore the metallicity and 
just focus on the age. To identify the final accurate age of 
each galaxy, we base our selections on the minimal \ 2 cr i~ 
terion. In particular, if all fitting ages for a single galaxy 
surpass the age of our universe, it is not convincing as 
there are several reasons causing this situation such as 
the low resolution of spectra and model dependence. Af- 
ter removing the result whose cij + cr, is greater than 16 
Gyr, where a, is the age and Oi is the age error, LRG 
ages are obtained and the age distribution is displayed 



Figure 3. Histogram of age spread information. This diagram 
shows the age result of 17832 LRGs with age interval of 1.0 Gyr. 
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Figure 4. Age-redshift envelope obtained from the galaxies stud- 
ied in this paper. Open circles correspond to LRGs from SDSS 
DR7. The red line shows the age of the universe t(z) for a ACDM 
model with Q m = 0.29 and Hq = 69 km s _1 Mpc - 1 while the blue 
line indicates t(z) + 3Gyr and is for reference only. A clear trend 
is present: galaxies age as redshift decreases. 

in Figure Figure 2] shows the ages of galaxies as a 
function of redshift. For clarity, the ages less than 7 Gyr 
have not been plotted. 

4. THE DETERMINATION OF OHD 

By using the age-redshift relation for SDSS LRGs, we 
now determine the Hubble parameter H(z). This proce- 
dure relies upon pinning down the 'edge' of the galaxy 
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Table 1 

Result of fitting 







The result of the first sample 




The result of the second sample 


n 


4 


5 


6 


7 


8 


5 


6 


7 


8 


9 


>x 


0.900i2 


0.94452 


0.99009 


0.98043 


0.99315 


0.96903 


0.91520 


0.99726 


0.99203 


0.97960 


SSE/n 


0.0526 


0.0758 


0.0492 


0.1064 


0.0947 


0.2506 


0.9640 


0.3193 


0.8018 


1.5750 


slope 3 


-9.53 


-f2.54 


-13.23 


-15.77 


-15.50 


-13.42 


-6.61 


-12.69 


-12.46 


-11.08 


°~slop 


6.07 


5.49 


3.76 


4.88 


4.33 


6.53 


9.62 


4.85 


6.74 


8.24 


—o slop /slope 


0.64 


0.44 


0.28 


0.31 


0.28 


0.49 


1.45 


0.38 


0.54 


0.74 


H b 


96.5 


72.8 


69.0 C 


57.9 


58.9 


64.8 


131.5 


68.6 C 


69.9 


78.8 






The 


result of the 


third sample 




The result 


of the forth sample 




n 


4 


5 


6 


7 


8 


10 


11 


12 


13 


14 




0.94508 


0.99f94 


0.99876 


0.99867 


0.99971 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


SSE/n 


0.0282 


0.0198 


0.0f69 


0.0338 


0.0309 


0.20844 


0.1130 


0.1752 


0.1631 


0.2109 


slope 3. 


-If. 40 


-14.44 


-11.14 


-11.83 


-12.12 


-11.70 


-7.20 


-8.61 


-7.34 


-7.81 


°~sLop 


12.33 


6.19 


4.53 


5.86 


5.34 


4.75 


3.09 


3.55 


3.15 


3.38 


—° slop/ 'slope 


1.08 


0.43 


0.41 


0.50 


0.44 


0.41 


0.43 


0.41 


0.43 


0.43 


H h 


71.4 


56.5 


72.9 C 


68.7 


67.0 


65.4 


106.3 


88.8° 


104.1 


98.0 



a slope — dt/dz, where t is the age of galaxies in the unit of Gyr. 

k H(z) is in units of km s -1 Mpc -1 . 

c Four new OHD obtained in this paper. 

age distribution in different redshift intervals. 

The main trunk of our produce is to fit the data t(zi), 
which is the oldest age at redshift Zj with straight line and 
its slope dz/dt is related to the Hubble parameter at an 
effective redshift by H(z e ff — —[1/(1 + z e f f)](dt/ 'dz)~ x 
where z e // = (z max + z min )/2. 

Mathematically, z rnax — z m i n should not be too great, 
so we divide these data points into four parts, each of 
which is called, expediently, a 'subsample' based on red- 
shift which are 0.03 < z < 0.11, 0.08 < z < 0.16, 0.16 < 
z < 0.24, and 0.24 < z < 0.32. Then we fit each 'sub- 
sample' by dividing it into suitably large redshift bins 
and get four new OHD finally. 

However, a problem is how to determine the number 
of bins, n. Practically, because of the incomplete LRGs 
sample, the age obtained is perhaps not the really oldest 
at that redshift point, that is to say, the age selected is 
perhaps the 'fake oldest'. We should try to avoid picking 
up 'fake oldest' age for the accuracy of our result, but 
too large n may increase the risk of picking up the 'fake 
oldest'. To indicate this risk, SSE/n (average of the 
Sum of Squares for Error), which is proportional to this 
risk, is calculated. Of course the number of 'n' should 
not yet be too small in order to limit the statistical er- 
ror which is proportional to the value of a s i ope / slope. 
So we divide each 'subsample' into 3, 4, • • • , 20 bins 
respectively, and every time we divide, we calculate the 
parameter A = (SSE/n, a s i ope / slope, 0P> X 2 ) • Finally, 
we compare all the parameters of A to get the value of 
'n'. But, there is something more to say, in the third 
and forth 'subsample', the error of the age is so big that 
the P> x ? is not such a good parameter to represent the 
goodness of fit. So in these 'subsamples', we consider 
the other two parameters more than the P> x ^ when we 
choose the result. Especially for the forth 'subsample', 
it is too large relative errors of slope that lead us to 
n > 10 (see in Tabled]). In addition, the optimal fitting 
of each part is also plotted in Figure [5] Finally we get 
four new OHD totally, which are -ff(z)=69.0± 19.6 km 

1 P> x ? is the probability of obtaining equal or greater value of 
the reduced x 2 if the *( 2 ) points are truly lying on a straight line 



Table 2 





H(z) measurements and their 


srrors 


z 


H(z) 
(km s _1 Mpc) 


o-h(z) 
(km s" 1 Mpc) 


Reference 


0.090 


69 


12 


1 


0.170 


83 


8 


1 


0.270 


77 


14 


1 


0.400 


95 


17 


1 


0.900 


117 


23 


1 


1.300 


168 


17 


1 


1.430 


177 


18 


1 


1.530 


110 


11 


1 


1.750 


202 


40 


1 


0.480 


97 


62 


2 


0.880 


90 


10 


2 


0.179 


75 


4 


3 


0.199 


75 


5 


3 


0.352 


83 


14 


3 


0.593 


101 


13 


3 


0.680 


92 


8 


3 


0.781 


105 


12 


3 


0.875 


125 


17 


3 


1.037 


154 


20 


3 


0.24 


79.69 


3.32 


i 


0.43 


86.45 


3.27 


4 


0.07 


69.0 


19. 6 a 




0.12 


68.6 


26. 2 a 




0.20 


72.9 


29. 6 a 




0.28 


88.8 


36.6 a 




References. — l[Simon ct ah] ]2005l) 2IStern et alj 



a New OHD points obtained in section HI 

s^Mpc" 1 at z = 0.07, (z)=68.6 ± 26.2 km s" 1 Mpc" 1 
at z = 0.12, H(z)=72.9 ± 29.6 km s^Mpc" 1 at z = 0.2 
and H(z)=88.8 ± 36.6 km s^Mpc" 1 at z = 0.28. Table 
[5] summarizes all the available OHD including the four 
new ones, and displayed in Figure H] 

5. COSMOLOGICAL CONSTRAINTS FROM OHD 

According to the Einstein gravity and Friedmann- 
Robertson- Walker (FRW) metric, H(z) read theoreti- 
cally: 



h{z) = H ox /n m (i + zf + n A + n k {i + z y 



(2) 



Hubble parameter 



5 




1 Fitting Result for the First Sample 
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Figure 5. The optimal fitted results. The solid line represents the best fitting for each sample. In the first sample 
Zmax = 0.109, and z e ff = 0.07. In the second sample, z m i n = 0.090, z max = 0.156, 2 e // = 0.12. In the third sample, 
z-max = 0.236, z e r t = 0.20, and in the forth sample, z m i n = 0.243, z ma x = 0.315, Z e ff = 0.28. 



0.033, 
0.170, 



where 
4.183 



Q r is ignored as it is too tinny at present (rs 
: 10~ 5 /i~ 2 ), and O m + Q\ + fi* = 1. Parame- 
ft\ and ilk parameterize the density of matter, 
cosmological constant and curvature density. The best 
fit parameters of the model via OHD are determined by 
minimizing the % 2 : 



N 



XoHu(Pmodel) - ^ 



(H obs ( Zi ) - H th {Zi)f 



(3) 



obs.i 



where p mo dei is a free parameter. Based on the 
basic equations of ACDM model, we have the fol- 
lowing two models: for a flat ACDM(57/ £ = 0), 



H th (z) = H oy /n m (l + z) 3 + (1 - n m ) with 
Pflat = (Cl m ,H ); for a non-flat ACDM, H th (z) = 



0(1 + ^) 2 , 
both models, 
as L oc exp 



With p„on-flat = (^A , O m , H ) . For 

the likelihood function can be written 

{-x 2 OHD m- 

We use MCMC method to calculate the likelihood in 
the independent parameter space. The Markov chains 
are gener ated and analyze d via the Python MCMC code 
- Pymc (|Patil et al.ll2010f) . The parameter Hq and the 
density parameter $7 = (Q m , ^a) are treatedas indepen- 
dent parameters in these two models (|Weill2010l ) , and all 
of the prior distributions of them are set as the uniform 
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Figure 6. All available OHD points. The dashed line plots the 
theoretical Hubble parameter Hfid as a function of z from the 
spatially flat ACDM model with Q m = 0.3, Qa = 0.7, and H = 
72 km s _1 Mpc _1 . 



Table 3 

The Prior of Model Parameters 



Model Parameters Prior Distrubution 



H 



Uniform(0.0,1.5) a 
Uniform(0.0,2.5) a 
Uniform(50,100) a 



a Uniform (lower limit, upper limit) stand for 
an uniform distribution in the interval [lower 
limit, upper limit], and the same below. 

distribution (Table [3]). 

Figure [7] and Figure [8] show the one-dimensional 
marginalized probability distribution for each parame- 
ter at the diagonal positions and the two-dimensional 
marginalized confidence regions at the non-diagonal po- 
sitions of the flat ACDM model and non-flat ACDM 
model respectively. The 3cr-s confidence regions 
(68.7%, 95.45%, 99.73% ) for each parameter are cal- 
culated detailedly. The best fit parameters are pointed 
out in the contour plots with a long vertical line and a 
long horizonal line across. 

To compare the constraining results of the 'new ver- 
sion OHD' within the new points we obtained and the 
'old version OHD', Figure[7]and Figure [8] show both two 
results of the same model with different line colors to 
distinguish: the red lines referring to the 'new version 
OHD' and the blue ones means the 'old version OHD'. 

The constraints in Figure [7] and Figure [5] show the well 
performances on constraining ACDM. On the flat-ACDM 
case, what we obtain of 'old version OHD' and 'new ver- 
sion OHD' are shown in Tabled] and for a non-flat sce- 
nario, the constraints are shown in Table [5] The new 
added-in points make amendments to the 'old version 
OHD'. In Fig [7] and Fig [5] the 1-D marginalized con- 
straints are more stringent than the old one and for the 
same model the contour plots against the same parame- 
ters pair have more shrunk confidence region at the same 
level. Since the number of new added-on points is so 
small, the degree of shrinking is weak as well. 

6. CONCLUSION AND DISCUSSION 



In this paper we present our measurement of four new 
OHD points from the ages of passively-evolving galax- 
ies at z < 0.4. A large sample of LRGs spectra has 
been dated by us and a age-redshift relation is obtained 
and displayed. By computing the relative ages of these 
LRGs, we give the final results on OHD. Combining them 
with other available 21 OHD points, we constrain cos- 
mological parameters using the update samples of OHD. 
The new added-in points data make amendments to the 
'old version OHD'. The 1-D marginalized constraints is 
more stringent than the old one and for the same model 
the contour plots against the same parameters pair have 
more shrunk confidence region at the same level. Since 
the number of new added-on points is so small, the de- 
gree of shrinking is weak as well. 

As our selected LRGs are incomplete and unable to 
cover all old-age galaxies in the universe, there would 
bring some difficulty in tracing the 'cosmic chronome- 
ters', that is to say, oldest age in our sample may not 
represent cosmic age at that redshift. It would be im- 
proved through future redshift surveys of observation and 
lead to successful realization of differential age method. 
In addition, as high SNR is essential for a precise fitting 
which determines the age of galaxies, we suggest that the 
accuracy of OHD would be improved if the SNR of spec- 
tra would increase. We employ ULySS to reconstruct the 
stellar population of galaxies. Though the robustness of 
ULySS has been illustrated, the fitting minimal may still 
be local because of the limitation of groups of initial val- 
ues. There has been significant advancement in modeling 
stellar populations of LRG galaxies and we expect these 
to improve the accuracy in this process in the future. 

The number of OHD is still scarce currently com- 
pared with SNIa data sets. OHD plays almost the same 
role as that of SNIa for the joint constraints on the 
ACDM model. The number of OHD points would ex- 
tended in further decades with more and deeper obser- 
vation of galaxies and at that time OHD set lone is 
capable to be u sed in place of current SNIa data sets 
(|Ma et al.1 120111 ). Fortunately, we have seen that BOSS 
project and proposed Sandage-Loeb observational plan 
(jCorasaniti et al.ll2007D can be used to extend our knowl- 
edge of cosmic expansion into the even deeper redshift. 
Finally, we think it is reasonable to expect that OHD 
will complement SNIa, BAO and weak lensing and help 
us getting more knowledge of the evolution history of our 
universe. 
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Figure 7. OHD Constraining results of flat ACDM. The diagonal plots show the 1-D PDF and the vertical lines from center to edge 
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levels, and the open blue circle point out the best fit point by 'old version OHD' while the red cross shows the ones via 'new version OHD'. 
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Table 4 

Marginalized Constraints on H m ,QA and Hq for flat ACDM Model 



Model Parameters 


Old Version OHD 




New Version OHD 
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+0.418484 
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n A 
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+2.154653 +4.570973 +6.899727 
-2.860217 -5.446933 -8.114114 



a The errors are la error, 2a error and 3a error, similarly hereinafter 
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Table 5 

Marginalized Constraints on Q m , Qjy, and Hq for non-flat ACDM Model 



Model Parameters 


Old Version OHD 




New Version OHD 
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+0.504268 
-0.304774 


302706+ 180116 
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-0.296992 
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-0.296992 


n A 
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+0.865635 
-0.702056 


0.658208t°; 3 3 6 5 6 4o5? 
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+0.904253 
-0.648651 


H ( ,(km s^Mpc- 1 ) 
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